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A new rescaling of the vorticity moments and their growth terms is used to characterise 
the evolution of anti-paraUel vortices governed by the 3D Euler equations. To suppress un- 
physical instabilities, the initial condition uses a balanced profile for the initial magnitude 
of vorticity along with a new algorithm for the initial vorticity direction. The new analysis 
uses a new adaptation to the Euler equations of a rescaling of the vorticity moments devel- 
oped for Navier-Stokes analysis. All rescaled moments grow in time, with the lower-order 
moments bounding the higher-order moments from above, consistent with new results from 
several Navier-Stokes calculations. Furthermore, if, as an inviscid flow evolves, this ordering 
is assumed to hold, then a singular upper bound on the growth of these moments can be 
used to provide a prediction of power law growth to compare against. There is a signifi- 
cant period where the growth of the highest moments converges to these singular bounds, 
demonstrating a tie between the strongest nonlinear growth and how the rescaled vorticity 
moments are ordered. The logarithmic growth of all the moments are calculated directly 
and the estimated singular times for the different Dm converge to a common value for the 
simulation in the best domain. 

PACS numbers: 

To appear in the Procedia lUTAM volume of papers for Topological Fluid Dynamics II under 
The growth of vorticity moments in the Euler equations 

I. BACKGROUND 

Two unresolved issues that have limited the application of numerics to the vortex dynamics and 
regularity questions of the three-dimensional Euler equations have been the inadequate analysis 
tools and the difficulties in specifying reproducible initial conditions. The existing analysis tools 
are unable to simultaneously cover the necessary range of scales in both space and time, while 
existing methods for mapping vortex tubes onto Eulerian meshes tend to generate ghost images 
unless ad hoc massaging is applied. This has led to weak and conflicting conclusions that depend 
upon the numerical method used and the choice of analysis that is applied to the results. 

To address these problems, this paper introduces an improved initialization for curved vortex 
tubes following an arbitrary trajectory and new analysis that is based upon higher-order vorticity 
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moments, and then appUes these to simulations of interacting anti-parallel vortices. The new 
initialisation suppresses core instabilities, which eliminates the ghost vortices found in earlier work 
|H] and discussed by [5]. Furthermore, the new trajectory algorithm allows the evolution of vortices 
with the same local perturbation, but different lengths, to be compared. One example calculation 
is given in Fig. 1. The new analysis allows one to compare all orders of the vorticity moments and 
generates new bounds against which to compare the results, which leads to more robust conclusions. 

The new analysis tool is an adaptation to the inviscid Euler equations of a new rescaling of the 
vorticity moments for the viscous Navier-Stokes equations. The rescaling uses a new frequency taJq, 
plus scaling powers Om, to convert the standard 17^; or L^"*, vorticity moments, into the following 
Dm moments [H [5] : 

Dm = (wq ^r^m)"'" where = ( / 1^1^""^^ J ^ wq = w„ = u / L?' and = 2m/ {Am 

(1) 

For the Navier-Stokes equations, vjq is based upon the viscosity v and the characteristic large length 
scale L of the turbulence and the am are designed such that neighbouring Dm (t) and Dm+i (t) terms 
can be compared directly using Navier-Stokes vorticity moment inequalities. This is adapted to 
the inviscid case below. The new rescaling makes comparisons between all the moments of the 
vorticity possible, both analytically and numerically. 

Historically, only the two limiting Dm have been used for addressing regularity questions. These 
are the global mean square vorticity or enstrophy, rescaled here into Di, and the point- wise max- 
imum of vorticity ||w||oo) rescaled into D^o- This is in part because analysis of the inequalities 
relating the intermediate moments had never been done. The known importance of Di is for 
addressing Navier-Stokes regularity (see references in [3]), while for the Euler equations, possible 
singularities are controlled by the time integral of ||w||oo- That is, for the Euler equations, if 

t 

\\uj\\oodT < oo for all time t > , (2) 

then the Euler equations are regular [1]. authors] Beale, J.T. 

The importance of the Dm between these limits is that by taking their ratios, new criteria for 
the regularity of the Navier-Stokes equations can be found [6j. To demonstrate the usefulness 
of the Navier-Stokes Dmit), Fig. 2(left) shows their evolution using data from a viscous, anti- 
parallel reconnection calculation using the initial condition described below. A full discussion 
of the trends, with lower order bounding higher order and convergence as m increases, is being 
prepared for publication. 
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To adapt this rescaling to vorticity moments of the inviscid Euler equations a non-viscous 
replacement for the scahng frequency vuq in ([T]) is needed. The inviscid modification chosen here 
defines wq using the circulation of the vortices T instead of the viscosity v. That is: zuq = wr = 

r/L2. 

For m < oo, a computational advantage of using these inviscid Dm in numerical analysis of the 
Euler equations is that they and their time derivatives dDm/dt can be determined at run-time and 
then compared as functions of time to integrals suggested by mathematical analysis. Furthermore, 
from the inverses of the logarithmic time derivatives {dlog{Dm) / dt)~^ = Dm/{dDm/dt), one can 
estimate the type of power-law singular growth using simple time differences [2] or, if it is assumed 
that the — )• a{Tm — t), running estimates of the Tm{t) can be made without using time 
differences. These running estimates will be used in the final test for singular growth using data 
from the best of the new anti-parallel Euler calculations. 

When these new Euler simulations were begun, the modest goal was to explain the type and 
strength of the convergence of the Dm moments in an early period of the Navier-Stokes calculation. 
The desired comparison period would be up to the beginning of the first vortex reconnection event 
at t = 16, shown in Fig. [Tj^right). Before t = 16, the viscous effects in the Navier-Stokes calculation 
should be negligible and the nonlinear terms, shared with the Euler equations, should dominate. 
The two frames in Fig. [2] are used to compare these Dm trends for the two Navier-Stokes and Euler 
calculations being highlighted here. 

The observed ordering of the Dm{t) for the Navier-Stokes simulations in Fig. ^Mt), plus 
the period of extended singular growth of the Dm^ Euler moments in Fig. [2| right) then led to 
Fig. |3j This figure addresses the question of whether the Euler equations have a singularity sub- 
ject] singularity using a new set of numerically determined time integrals and analysis of logarithmic 
time derivatives found at run-time. The new time integrals come from mathematical analysis of 
the Euler equations that assumes a priori that the order of the Dm seen numerically will hold for 
all time, as indicated by Fig. [2]^right). 

All the calculations are, fundamentally, in periodic computational domains, with symmetries 
used to decrease the data and time needed to do the calculations. Several filtered/dealiased pseu- 
dospectral methods have been tested and described previously [2j. The method chosen for the 
calculations here is a combination of the 2 /3rds dealiasing rule plus a 36th order filtering method 
that was first introduced without dealiasing [8j. authorsjKerr, R.M. authors]Bustamante, M.D. 
authors]Hou, T.Y. 

The principle axes are: x is the direction of propagation of the vortex pair, y is in the primary 
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direction of the vortices, and z is the direction between the vortices. The computed domain size 
is Lx >i Ly X Lz, while the fuhy-periodic domain would be in Lx x 2Ly x 2Lz- Domain sizes and 
meshes are given in Table 1. Referring to the initial condition in Fig. [T| the y = symmetry plane 
with the maximum perturbation will be called the perturbation plane and the z = symmetry 
plane between the vortices will be called the dividing plane. 

The paper is organised as follows. First, the new initial condition is described briefly. Next, the 
re-scaling of the vorticity moments for the Navier-Stokes and Euler equations is discussed further 
and applied to the new calculations, from which a new ordering for the Dm{t) moments is found. 
This ordering has been found for all times for both the viscous and inviscid cases. After the new 
ordering is established, new upper bounds on the growth of the moments in the Euler equations 
are found and applied to the inviscid Euler solutions. Next, the logarithmic time derivatives of the 
from the Euler calculation are used to give running estimates of the singular times, labeled 
Tm{t)- It is found that for m > 1, these estimated times converge. That is, all the Tm{t) — )• Tc- 
Finally, there is some discussion of additional diagnostics that are now being collected, such as 
the curvature of the vortex lines, that will be needed if we are going to understand why the Euler 
equations can obey singular scaling laws for extended periods. 

II. INITIAL CONDITION 

At meeting on the Euler equations in 2007 in Aussois, France, one topic was results from 
direct numerical simulations that addressed the question of regularity of the Euler equations. 
The conclusions of the two anti-parallel calculations [21 E] were different, even though both were 
nominally using initial conditions similar to [llj. authors]Kerr, R.M. authors] Bustamante, M.D. 
authors]Hou, T.Y. Clearly, the prescription in [11] was flawed. These flaws have now been identified 
and will be described in detail in another paper. The three primary elements of the new initial 
condition are these: 

• A new profile of the vorticity distribution in the core that is based upon the Rosenhead 
regularisation of a 2D point vortex and is similar to the two-dimensional density profiles 
used for quantum Gross-Pitaevskii calculations |13j . 

• A new direction algorithm that, for a given {xi,yj, Zk) on the three-dimensional grid, begins 
by finding the nearest position {xs,ys-,Zs) on the given analytic trajectory. The distance 
used in the profile function for finding \uj\{xi^yj^ Zk) is r = \{xi^yj, z^) — {xs^ys-, Zs)\ and 
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the direction of the vorticity at the points {xi,yj,Zk) is given by the tangent of the chosen 
trajectory at {xs,ys,Zs). 

• Making the vortices very, very long to minimise boundary effects. 

The resulting profile has been used for anti-parallel vertical vortices in a stratified fluid, the 
anti-parallel unstratified Navier-Stokes calculations mentioned here, and now anti-parallel Euler 
vortices. In each case, unphysical initial instabilities due to small-scale inbalances have been sup- 
pressed, a cleaner and stronger larger-scale instability has been identified, and, where appropriate, 
a transition to sustained turbulence forms from the vortex interactions where none had been seen in 
earlier work. Unlike in earlier work [21 [TT] , no extra massaging or squeezing of the initial condition 
is needed to ensure that there is only one sign of the vorticity in the calculated y = perturbation 
plane. 

The computational procedure is as follows: First, the vortex is initialised on a modest mesh, 
which is then put onto a much larger computational mesh by adding zeros at the higher wavenum- 
bers. The calculation then proceeds on this mesh until, by comparing results on different meshes, 
the collapse has progressed to the point where the calculation would soon be underresolved. Then 
the calculation is remeshed onto a finer mesh. Two remeshings are typically needed to reach the 
final mesh at the final times. The calculations used for the current study are given in the Table 1. 

The initial and evolved Navier-Stokes vorticity isosurfaces in Fig. [l] apply to both the Navier- 
Stokes and Euler calculations because viscous dissipation for the Navier-Stokes case at t = 16 has 
been minimal. The insets show the upper/left quarter domain near the y = perturbation plane, 
with the t = inset showing that the initial vortex tube has a circular cross-section of constant 
width along its entire length. 

These figures can be compared with similar stages in the evolution of anti-parallel quantum 
vortices in [13] and to the cover illustration in pL2j , which shows how the vortex lines twist as they 
extend from the y = perturbation plane. The vortices do more than twist. They actually bend 
back upon themselves until, near y = ±7, they are closer than the unperturbed original vortices for 
y > ±8. The possible significance of this bend and its curvatue will be discussed in the summary. 

III. NAVIER-STOKES INTERMITTENCY AND THE RESCALING VORTICITY 

MOMENTS 

subject] vorticity moments 
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T=0 v^O.001 1(01=0.38 maxl(ol=1.0032 T=16 v=0.001 lojl=0.77 maxlcol=8.7 




FIG. 1: Left: Very long, anti-parallel initial condition at t = 0. Right: t = 16 (Navier-Stokes). Insets 
show z > for < y < 10. 

A neglected topic in studying the Navier-Stokes equations is temporal intermittency, periods 
of intense activity, interspersed by relatively quiescent periods. One approach to characterising 
this type of intermittency is through higher-order strain ([S] = Sij = 0.5{dui/dxj + duj/dxi)) 
and vorticity (a; = V x u) moments, plus experimentally measurable single-point derivatives |17] . 
Numerically, convergent statistics for S^™" and u"^"^ with m=2 and 3 were obtained as early as 1985 

m- 

However, having only orders m=2 and 3 is insufficient for making theoretical comparisons and 
it has been impossible to get convergent statistics for the next higher-order moments for even the 
largest forced simulations [9j. The problem is two-fold. First, the difference between the higher- 
moments moments in the quiet periods and the intense periods can be huge, and second, these 
occur on the time-scale of the large-scale forcing for simulations that, due to their size, can only 
be run for a few of these characteristic timescales. 

Recently, Yeung, Donzis & Sreenivasan |18j authors] Yeung, P.K. have found that convergent 
statistics for their forced simulations can be obtained by taking ratios of the higher-order moments. 
While simultaneously, new mathematics has concluded that these ratios, rescaled in a manner 
consistent with inequalities for the time derivatives of the higher-order vorticity moments Dm 
[1], can give new insight into the Navier-Stokes singularity question 0, authors] Gibbon, J.D. as 
summarized in [14j . 

One analytic approach to answering whether the Navier-Stokes equations are regular or not 
starts by assuming that there are quiescent and intensely intermittent periods, called good and 
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FIG. 2: Left: Re = T/u = 4000; Dm Q fr om an anti-parallel calculation. Dm are ordered with lower-order 
bounding higher-order for all times. Especially note the the periods of steepest growth, t < 16 and t ~ 90 when the 
nonlinear terms dominate. Right: The inverses: D^{t), from the Euler calculation with a similar initial condition 
for the first period of sharp growth {t < 15). The hierarchy of D^{t) includes ti7/||a;||oo. The Dm^it) from the 
Navier-Stokes calculation are similar, but with a greater deviation from a linear form as time increases. 



bad, or possibly neutral [3]. Because the Dm{t) moments can be compared directly using vorticity 
moment inequalities, the new mathematics [5] is able to derive new bounds on the periods of the 
maximum growth of the Dm that can be compared to numerical results. 

Fig. [2]^left) shows how the are ordered for the new reconnection calculation, with the lower 
order Dm bounding the higher order D^. for all times: Dm+i{t) < Dm{t). This means that the 
different Dm never cross one-another and the definitions of bad and good periods are the same 
for all of the Dm- This ordering of lower-order above higher-order Dm was unexpected because 
it is opposite of how the original Om, without rescaling, are required to be ordered using Holder 
inequalities and is opposite to what would easily ensure regularity of the Navier-Stokes equations 
using the new bounds of [5]. authors] Gibbon, J.D. This ordering has now been identified in every 
Navier-Stokes simulation it has been tested against. A joint paper is in preparation and mentioned 
in 0. 

What governs the dynamics during these intense/bad and quiet/good periods? Fig. [T]^ right) 
shows the structures at the end of the most extreme growth of the Dm{t), up to t = 16, when the 
first Navier-Stokes vortex reconnection is forming. In a new Navier-Stokes reconnection paper it 
will be shown that all of the subsequent periods of intense growth of the higher-order Dm can be 
tied to how vortices are attracted and stretched just before reconnection events. 
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IV. RESCALED VORTICITY MOMENTS AND THE EULER EQUATIONS 

The strongest growth of the Navier-Stokes Dm in Fig. [2]^left) is before the first reconnection at 
t w 16, when viscous effects are neghgible and the nonhnear Euler dynamics are strongest. During 
this period, the growth of the normahsed enstrophy production, a skewness factor, is is up to three 
times the values typically determined in large Reynolds number experiments and simulations. In 
order to understand the origins of this period of growth, a new series of simulations of the Euler 
equations were begun that cover part of this period {t < 14.25). The particular calculation shown is 
just one of a series of new anti-parallel inviscid Euler calculations outlined in Table 1, all using the 
new profile and trajectory algorithms. The objective was to find, and confirm, whether a domain 
could be identified where the boundaries were not suppressing any growth of the maximum of 
vorticity || cj||qq. And then, determine if this calculation indicates singular behaviour, or not. 

It was found that changing the length of the domain in the y-direction had the greatest effect 
upon the growth of ||w||oo, with its growth being suppressed until Ly = Bvr was reached, the case 
labeled vllbzz. The vllg case with Ly = IGtt gave identical results to the vllbzz case. 

The choice of vjq = wy for the Euler Dm analysis was in part inspired by how the growth of 
enstrophy in an earlier Euler calulation [15J could be explained empirically by replacing viscosity 
u with the circulation T in the well-known inequality for the upper bound on enstrophy growth in 
Navier-Stokes:: {d/dt)^l < Ci{9\/ {v / L'^)f . This empirical guess can now be replaced by robust 
Euler bounds that use the Dm with wq = wy- 

Following the proof of Proposition 1 [6j, one starts with: 

2mL3Q^— < 2mL^ci,mKX\^Z (3) 



which, with some rearranging, becomes 



m+l 



,,«„<c.„, ^1 ni.. (4) 



Finally, upon substituting the definition of the Dm and pulling the out, one gets 
dt 



d ( Dm 1 \ ^"^ 

^-Dm < C2,m^o ( \ DI^ where Cm = \{^'m + 1) . (5) 



Once ujo = is chosen, then the inviscid Dm{t) can be compared for the Euler calculation. 
This has been done in same manner as in Fig. [2]^left) and shows the same ordering as in the 
Navier-Stokes case[Tn]. However, in order to include the m — t- oo limit, a better choice is to plot 
Dm^{t), which for m = oo gives D^ = zur/\\u)\\oci where \\uj\\oo is the sup(|a;|). A simple test for 
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singular behaviour is to compare 1/|| 

'^lloo against the power law consistent with the lower bound, 
for singular growth of ||u;||oo allowed by Q. That is: 

\\uj\\oo ^ (Too - or Dq^^ = -ajr/llwlloo ~ tur(Too - t) . (6) 

Under this test, the sign of singularity growth would be finding that — t- linearly. 

The D^{t) are plotted in Fig. [2|right). Compared in this way, as m becomes large, the 
nearly match uJriToo — t) as t increases. However, the growth of || 

^||oo appears to tail off of this 

behaviour as t — )• 15. So, to claim singular growth, another independently calculated diagnostic is 
needed to confirm the trends seen in Fig. [2](right). 

If the only diagnostic for singular growth is ||w||oo, then an appropriate secondary diagnostic 
could be a = dlog HwHoo/c^i, the logarithmic time derivative of ||ct;||oo- In principle one should 
determine a from the vortex stretching exactly at the position of || However, to get the 

stretching at the exact position of ||w||oo, which lies between the mesh-points in physical space, 
requires interpolation, which can be both difficult and inaccurate. In practice, the only stretching 
diagnostic that did not have grid-induced oscillations and was near, but not at, the position of 
1 1 a; I loo, was to take the maximum of the vortex stretching on the perturbation plane [llj . 

Using the Dm resolves this problem. The trick is to calculate both the Dm and their time 
derivatives {d/dt)Dm at run-time, a simple matter of programming compared to the interpolations 
needed for determining a at ||w||oo- Furthermore, Fig. [2]^ right) shows that as time and m increase, 
the D^{t) — )• D^{t) = wr/\\uj\\oo ■ Therefore, for large m, the secondary diagnostics equivalent to 
a are the logarithmic time derivatives of the Dm (t) , which are used below to define the estimated 
singular times Tm{t) 



Numerical analysis using new Euler integrals 

With the added assumption that the Dm+i/Dm are always bounded, as demonstrated by Fig. 
[2|^ right), let us use the bound in ([s]) to help us write new Euler bounds that can be tested numer- 
ically. 

For general m, let us begin by rewriting ([s]) as 

-^-Dm^ < C^ror ^^^^^^ ' ^Si'^g Fm{t) = Cm J "^V (^ ^^^ ^ ^'^ One gCtS D:^^ < C2,mFm{t) 

This focuses our attention upon the integrals on the right-hand-side, which are plotted in Fig. 
(sfleft). In this fi gure, the upper bound, based upon the integral of D^ , grows linearly. If obeyed 
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exactly, this would imply that the solutions are singular. However, since this is only an upper 
bound, another test is needed. 

This final test will be a diagnostic coming from the logarithmic time derivatives of the D^. 
Two assumptions are made. First, that D'^{t) ~ (Tm — and second an assumption on the 

7m. Applying {dlog D'^/dt)^^ to the first assumption, one gets 

idlogDl/dt)-'=^-\Tm-t). (8) 

By applying this assumption to time differences of the {dlog D'^/dt)~^ , one could get running 
estimates of both the and the T^. However, since as m increases the curves generated by the 
D'^(t) are becoming linear in both Fig. [2|right) and Fig. [sj^left), that is — ^ !> the best way to 
find running estimates of the Tmif) is to make 7m = 1 an added assumption and use 

Tra{t) = {d\ogDl,/dt)-^ +t. (9) 

The result is in Fig. |3||right). For t > 12, the estimated m > 3 singular times Tm{t) are 
beginning to converge. This is shown more clearly by adding a. t = t curve and extending the 
computed Tm with linear extensions based on the Tm at the last two times computed. If there is 
a singularity of the Euler equations for this initial condition, then they should all cross the t = t 
line at the same time. Which they do. 



TABLE I: Domains and sequences of meshes used to resolve. 



Domain 


label 


Mesh 1 and Mesh 3 


At 




Mesh 2 and Mesh 4 




At 




Stt X Stt X 2n 


Vila 


512 X 256 X 512 


t = 0- 


- 12 


512 X 256 X 1024 


t 


= 8- 


- 15 


StT X StT X 277 


vlla 


1024 X 512 X 4096 


< = 10 - 


13.25 


1024 X 512 X 2048 


t 


= i- 


- 15 


4tt X in X 277 


v71 


512 X 512 X 1024 


t = 0- 


- 12 


1024 X 512 X 2048 


t = 


- 12 - 


13.5 


4tt X in X 2n 


v71 


1024 X 512 X 4096 


t = 13.5 - 


- 14.25 










in X 6n X in 


vllby 


512 X 256 X 2048 


t = - 


- 12 


1024 X 512 X 2048 


t 


= 8- 


- 14 


in X in X in 


vllbx 


1024 X 512 X 2048 


t = Q- 


- 14 


1024 X 512 X 4096 


t = 


11 - 


14.25 


in X 8n X in 


vllbz 


512 X 512 X 1024 


t = 0- 


- 12 


1024 X 1024 X 4096 


t = 


12- 


13.75 


in X 8n X 2n 


vllbzz 


512 X 512 X 1024 


t = 0- 


- 12 


1024 X 512 X 2048 


t = 


: 12 - 


13.5 


in X 8n X 2n 


vllbzz 


1024 X 512 X 4096 


t = 13.5 - 


- 14.25 










in X 16n X 2n 


vllg 


512 X 1024 X 1024 


t = Q- 


- 12 


512 X 1024 X 2048 


t = 


12 - 


14.25 


in X 16n X 2n 


vllg 


1024 X 2048 X 2048 


i = 10 - 


- 14 


1024 X 2048 X 4096 


t = 


: 13- 


14.5 
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FIG. 3: Left: J ^^{Djn^i/ D^)^'^ dt ([7]) for an anti-parallel Eulcr calculation. The observed ordering and 
linear increase as m oo would permit at least super-exponential growth of the D„i for both Euler and 
Navier-Stokes up to t = 14.5. Right: Estimated singular time from different T„i{t) = {d\ogDj^/dt)^^ + 
t. Only m odd are shown to reduce clutter and a curve with t = t in added to clarify where the Tm(t) are 
heading. Linear extrapolations to t — 15.75 of the m > 1 curves, based on the last two values, are shown 
with the dashed lines. As well as could be expected, these extrapolations all appear to be crossing the 
t — t line at about Tc ~ 15.8. For t > 12 (and excluding m = 1), the Tm{t) are ordered. Going from 
underestimating the {m — 3, 5) to overestimating T^. (m = 7, 9). 



V. SUMMARY 



A new approach to rescaling vorticity moments, the Dm, has been used for the analysis of 
new Navier-Stokes and Euler calculations. The Dm have the following favourable analytic and 
numerical properties: In mathematical analysis, neighbouring orders can be compared using their 
time derivatives inequalities [6j. In numerical analysis, their values, time derivatives and thus their 
logarithmic times derivatives can be determined continuously and compared. 

The numerical comparisons have revealed an unexpected hierarchy where the lower-order Dm 
bound the higher-order Dm, for all times and for both sets of calculations. This ordering of the 
Dm was unexpected for two reasons. First, it is opposite to the required Holder ordering of the 
and second, it is opposite to an ordering that would immediately imply that the Navier-Stokes 
equations are regular for all times. 

Furthermore, the period of strongest growth and alignment of the Dm{t) occurs when the 
normalized enstrophy growth for the Navier-Stokes calculations is strongest. These observations 
led to the secondary goal of the Euler calculations. This is to use the Dm, and their logarithmic 
growth rates, to determine whether these Euler calculations are consistent with the formation a 
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finite-time singularity. 

The new Euler calculations are consistent with the formation a finite-time singularity in the 
sense that each of the higher-order moments in Fig. [2] show singular trends for a longer period than 
any previous Euler calculation. This includes D'^, the rescaled singular maximum of the vorticity 
1 1 w I loo- Critical to achieving this extended period of singular Euler growth is using vortices that are 
not subject to internal instabilities and domains that are longer than in any earlier work. However, 
using the new profile and direction algorithm is not enough. For the smaller domains listed in 
Table 1, the singular growth saturates early, as in some earlier work |8j. The importance of the 
longer domain is that it allows the full effect of the curvature and torsion of the vortices to manifest 
itself. It was not until the length of the domains was Ly = Svr that the extended period of singular 
growth appeared. 

With the new data set and new results, a number of outstanding questions will soon be ad- 
dressed. One is determining the role of the curvature of the vortex lines. New analysis shows that 
the curvature of the vortex lines near HwHqq is small, and therefore contributes little to the local 
y = vortex stretching. This suggests that the stretching is coming from the strong curvature and 
looping seen in Fig. [l] (right) at t = 16 for y = ±5. Even though this structure is forming far from 
the position of ||w||oo on the y = perturbation plane, its effects upon the growth on that plane 
are surprisingly strong. 

This success with finding the curvature suggests that are further conditions on properties derived 
from the direction of the vorticity can be determined and tested against the growth of ||u;||oo .The 
new data should be capable of testing these proposed constraints and identifying how the position 
of the maximum of vorticity moves with respect to the Lagrangian flow. This relative motion is 
non-zero [? ] and could be important for determining how the circulation in the perturbation 
plane becomes divided into a head and flattened tail, where the head is the part of the circulation 
that reconnects if there is viscosity. 
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